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ABSTRACT 

The Optica) Properties of Astronomical Silicates with Infrared Techniques program utilizes multiple instruments 
to provide spectral data over a wide range of temperatures and wavelengths. Experimental methods include Vector 
Network Analyzer and Fourier transform spectroscopy transmission, and reflection/scattering measurements. From 
this data, we can determine the optical parameters for the index of refraction, n, and the absorption coefficient, 
k. The analysis of the laboratory transmittance data for each sample type is based upon different mathematical 
models, which are applied to each data set according to their degree of coherence. Presented here are results from 
iron silicate dust grain analogs, in several sample preparations and at temperatures ranging from 5 to 300 K, across 
the infrared and millimeter portion of the spectrum (from 2.5 to 10,000 pm or 4000 to 1 cm -1 ). 

Key words: dust, extinction - infrared: ISM - methods: laboratory - submillimeter: ISM - techniques: spectroscopic 

Online-only material: color figures 


1. INTRODUCTION 

Interstellar dust is found in virtually every astrophysicai 
environment, including the solar system, star-forming regions, 
stellar debris disks, and distant galaxies. This dust plays an 
important role in the chemistry and radiative properties in 
these environments, and hence helps determine the evolution 
and ultimate state of these objects. It acts as a catalyst for 
reactions of molecular hydrogen, hydrocarbons, and water. Dust 
ejected from supemovae and stellar outflows of evolved stars can 
enter the interstellar medium (ISM), ultimately coalescing into 
the disks surrounding newly forming stars and planetesimals. 
Furthermore, dust has profound effects on the light traversing 
such regions. For example, dust concentrated in the plane of 
the galaxy is an effective barrier to visible and ultraviolet 
wavelengths, and dust grains surrounding stellar regions or 
in the ISM can redden higher frequencies of observed light 
by absorption, scattering, and long-wavelength re-emission of 
absorbed light. At least 30% or more of the energy emitted as 
starlight in the universe is re-radiated by the dust in the infrared 
(IR; Bernstein et al. 2002). 

There is little doubt that silicate materials contribute a 
substantial fraction of the total mass of interstellar dust. The 
spectral feature at 10/im was first observed by Gillett et al. 
(1968), and shortly afterward was identified as resulting from 
the presence of silicate grains by Woolf & New (1969). Since 
then, astronomers have attempted to use observations of the 
10 and 18/tm bands to determine which silicate species are 
responsible for the observed emission (Henning 2010a). 

While laboratory spectra of “amorphous” silicate materials 
currently available in the astrophysicai literature are typically 
those of glassy solids (e.g., Boudet et al. 2005; Henning & 
Mutschke 1997), astronomers have been lacking high quality 
laboratory data of analogs to astronomical dust grains (see Speck 
et al. 2011, Table 1 for a recent review of types of laboratory 
samples). Laboratory spectra for chaotic materials, such as the 
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ones developed by Nuth & Hecht ( 1 990), have been investigated 
over a limited spectral range (e.g., mid-IR wavelengths) at 
room temperature (Hallenbeck et al. 2000); however, it is 
essential that the data set for such analogs be extended to longer 
wavelengths of the IR spectrum and to cover a significant range 
of temperatures. The optical properties of silicate dust grains 
at low temperatures have been shown to differ from those at 
room temperatures (Day 1976; Zaikowsky et al. 1975) and 
evidence suggests possibly significant temperature-dependent 
variations of spectral feature shapes and positions (Chihara 
et al. 2001). New condensate spectral data collected over an 
extensive spectral range are important for the interpretation of 
existing observational data from missions such as Spitzer and 
Herschel . and will be critical for understanding observations 
from facilities such as SOFIA, James Webb Space Telescope 
C JWST ), and ALMA. 

The Optical Properties of Astronomical Silicates with In- 
frared Techniques (OPASI-T) program in progress at Goddard 
Space Right Center was designed to provide such data (Kinzer 
et al. 2010; Rinehart et al. 2011). The objective of OPASI-T is to 
measure the spectra of disordered silicates in order to constrain 
the complex optical parameters, the real part of the complex 
refractive index (n, referred to as the index of refraction) and 
the imaginary part (k, referred to as the absorption coefficient), 
which may vary as functions of dust temperature and crys- 
tallinity. Spectra are collected from the millimeter and far-IR 
wavelengths through the mid-IR, in order to connect the new 
data with existing laboratory and astronomical measurements in 
the mid-IR. 

2. EXPERIMENTAL METHODS 

Due to the frequency dependence of the optical parameters, 
the OPASI-T program makes use of different sample preparation 
methods, as well as multiple instruments, to provide laboratory 
data of iron silicate dust grain analogs. An extensive discussion 
of these methods and apparatus is presented in Rinehart et al. 
(20 1 1 ). which also included results on SiO samples used to 
validate our experimental techniques. 
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** Co 1 ndensadon "" A W~ (CFA) used .o manufacture nanometer-scale smoke panicles via die combustion of hydrogen gas 
that contain small amounts of silane and iron pentacarbonyl. Note that this figure was previously used and described in Nuth et al. (2002). 


2.7, Sample Production of Dust Analogs 

Previous and ongoing laboratory investigations of silicates 
have tried to fit spectra of terrestrial minerals, such as forsterite 
and enstatite, to measured astronomical emission spectra. Such 
materials are relatively easy to produce in large quantities and 
are reasonable first approximations to dust in astrophysical en- 
vironments (for reviews, see Henning 2010b; Speck et al. 2011). 
The crystalline and glassy minerals common in terrestrial en- 
vironments result from high temperatures, an oxidizing atmo- 
sphere, and the cooling process (slow cooling present in mag- 
matic systems produces ordered ciystals while rapid cooling 
produces well-ordered glasses). Silicate tetrahedra are the basic 
building blocks of all such materials. The tetrahedra of crys- 
talline structures form well-ordered planes with cations (e.g., 
Fe, Mg) filling specific vacancies in the matrix (and providing 
electrostatic neutrality), while glasses (low temperature) and 
melts (high temperature) have quasi-random orientations of the 
silicate tetrahedra. High doses of ionizing radiation can create 
defects within the crystal structure. High concentrations of such 
defects can have a considerable effect on the spectra (Day 1977; 
Demyk et al. 2001). 

Glassy materials can have ordered structure on microscopic 
scales, with a high level of disorder at macroscopic scales; as- 
tronomical dust grain, however, are largely disordered on all 
scales. Typically, silicon in a high temperature neutral or ox- 
idizing gas is in the form of SiO while metals are present 
as atoms (Fe), monoxides (A10), hydroxides (AIOH), or hy- 
drides (A1H). Solids form at veiy low pressures and kinetic pro- 
cesses “freeze out” solids before equilibrium is achieved. This 
produces solids with chaotic fundamental chemical structures 


that lack the simple silicate tetrahedral building blocks. These 
solids do not even achieve uniformity of the chemical ioniza- 
tion states of the silicon atom, let alone uniformity of cations 
in the solid. Although one may think of these solids as highly 
defect-rich, this can lead to the misinterpretation that the mate- 
rials might once have been highly ordered. In reality, if the gas 
cooled rapidly enough, then the solids are a close representation 
of a frozen gaseous state. These condensates may contain neutral 
atoms, and can easily adsorb SiO and H 2 O onto their surfaces; 
as the grains grow through additional adsorption, layers of these 
materials are “buried” within the grain. While other groups have 
also have also produced materials that show thermal annealing 
to an originally amorphous condensate (Speck et al. 2008), we 
are providing these data to test the hypothesis that grains ini- 
tially form as amorphous solids followed quickly by some level 
of thermal annealing. The more effective the annealing (higher 
temperature, longer times), the more crystalline the final grain 
population will be (Nuth et al. 2002). 

We produce astrophysically relevant silicate grain analogs 
within a dust generator constructed at Goddard Space Flight 
Center. The Condensation Flow Apparatus (CFA), shown in 
Figure 1, was designed to condense amorphous grains with 
chemical compositions determined stochastically by the com- 
position of the vapor from which they formed. It produces gram- 
quantities of aggregate samples of amorphous grains (Nuth 
et al. 2000a) by non-equilibrium gas-to-solid condensation in 
M-SiO-H 2 -02 vapors (M: Mg, Fe, Al, Ca, or combinations 
thereof; Nuth et al. 2000b, 2002). The atmosphere in the appara- 
tus has a total pressure of 90 Torr and is dominated by hydrogen 
at temperatures in the 500-1500 K range. Condensable species 
form very rapidly via combustion from gas-phase precursors 
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Tbblel 

Sample Types and Wavenumber Coverage 


Sample Type 

Spectral Types 

Wavenumber Range 

Temperature 

Packed dust in waveguide holder 

Transmission 

— 1 citT 1 

295 K 

Packed dust in aluminum holders 

Transmission 

30-250 cm" 1 

5-295 K 

Dust in PH matrix 

Transmission 

40-650 cm" 1 

5-295 K 

Dust in KBr matrix 

Transmission 

400-4000 cm" 1 

295 K 

Dust+acctone on aluminum disk 

Reflection 

100-650 cm" 1 

5-295 K 

Dust+acetone on aluminum disk 

Reflection 

400-4000 cm" 1 

295 K 


such as silane (SiH 4 ) and iron pentacarbonyl [Fe(CO) 5 ]. These 
species typically constitute <10% of the total gas input to the 
system. The oxidant, usually pure oxygen (although nitrous ox- 
ide has also been used), is introduced just before the furnace. 
Finally, a volatile metal, such as Mg or Ca, to be combined with 
the silicate, can be placed into the furnace within a graphite cru- 
cible. The furnace temperature controls the vapor pressure of the 
metal. The flow velocity of the metal vapor through the furnace 
is typically on the order of 10-20 cm s~\ A grain formed at the 
high-temperature flame front near the furnace entrance spends 
less than a second within the furnace following nucleation and 
growth. The hot gas and fresh grains are rapidly quenched as 
they flow into a larger stainless steel chamber lined with an 
aluminum substrate at a temperature of 300-350 K. 

Condensation and growth within the dust generator is a 
stochastic, kinetically controlled process, often producing grains 
that are fluffy, open aggregates. Typical grains are on the order of 
—2-30 nm in radius and aggregates (clumps) frequently consist 
of hundreds to thousands of individual grains, each connected 
to only two or (rarely) three neighboring particles (Nuth et ah 
2000b). The dust generator was designed to condense amor- 
phous grains with chemical compositions determined stochasti- 
cally by the composition of the vapor from which they formed. 
However, the grain compositions match metastable eutectics in 
the appropriate phase diagrams. A surprising result of these ex- 
periments has been to show that iron silicates and magnesium 
silicates condense as separate grain populations even from a 
very well mixed vapor (Rietmeijer et ah 1999a). By measuring 
the properties of these simple grain compositions, we can ob- 
tain grain spectra representative of dust in astrophysical systems 
(Rietmeijer et ah 1999b). The size, morphology, and crystallo- 
graphic properties (i.e., amorphous or crystalline) of individ- 
ual grains were characterized by 200 keV analytical and high- 
resolution transmission electron microscopy (TEM); for details, 
see Rietmeijer et ah (2006). These TEM analyses showed that 
(1) the exact bulk composition of the condensing gas is un- 
known, though it is very close to the bulk composition of all 
solid grains in the resulting smoke (Nuth et ah 2000b, 2002), 
and (2) formation of hydrated condensate grains by transient 
water in the condensation chamber is extremely rare Rietmeijer 
et ah (2004). Note that due to the grain size being well below the 
Rayleigh limit, the grain shape should only affect the extinction, 
and not our measurements. 

In order to characterize the transmission and reflection of 
the material throughout the wavelength range covered in this 
program, several sample preparation techniques were used and 
are shown in Table 1 . For millimeter and far-IR transmission 
measurements, the iron silicate dust is packed into holders (e.g., 
the waveguide or aluminum holders for Fourier transform spec- 
troscopy, FTS). For reflection measurements, the iron silicate 
dust is diluted into acetone, then deposited onto an A1 disk 
for measurements. Transmission measurements in the mid-IR 
required the iron silicate dust to be embedded as an inclusion 


within a matrix material (e.g., KBrorPE). Matrix measurements 
can present challenges to modeling and interpretation, but sim- 
ilar challenges exist for all sample preparations, and embedded 
samples are best-suited to these materials at these wavelengths. 
The validity of the experimental technique was demonstrated 
using SiO, as discussed in Rinehart et al. (201 1). Each sample 
preparation technique is described in the measurement sections 
below. 

2.2. Waveguide Measurements 

At microwave frequencies, the physical size of the dust parti- 
cles is very small in comparison to the wavelength and the sam- 
ple appears homogenous to the incident radiation field. The ma- 
terial properties can be readily treated by the Maxwell-Gamett 
effective medium approximation, where the two phases (sil- 
icate and air) are randomly dispersed and form a separated 
grain composite mixture (Niklasson et al. 1 98 1 ). A Fabry-Perot 
waveguide resonator technique (Wollack et al. 2010) is used to 
characterize the material’s permeability and permittivity func- 
tions, as a function of frequency, with an Agilent E8364A PNA 
Series Vector Network Analyzer ( VNA). In this cavity measure- 
ment approach, the dust sample is packed in a waveguide cavity 
that is closed on each end with a thin window (i.e., — 1 2 /zm Kap- 
ton film, American Durafilm). The waveguide used here is an 
8^00 = 203.2 mm long section of WR28.0. The guide broadwall 
and height are 0?280 (7.112 mm) and tf.'140 (3.556 mm). The 
sample holder has standard UG-599 square waveguide cover 
flanges which are 01115 (19.05 mm) on a side. This smooth in- 
terface is used to secure the thin window material against the 
waveguide flange that defines the VNA’s calibration reference 
plane. To enable repeatable and precise dielectric parameter ex- 
traction from the observations, the sample packed into the holder 
must have a high degree of homogeneity. The desired unifor- 
mity was achieved by lightly tapping and adding material until 
the sample holder is filled; this procedure helps to eliminate 
voids and remove trapped air, which can lead to complications 
in interpretation of the observed spectra. 

2.5. FTS Transmission Measurements 

Samples for FTS transmission measurements are prepared in 
two ways. For longer wavelengths (far-IR), powders are packed 
within a custom-made aluminum holding cavity (10 mm diam- 
eters with varying cavity lengths) equipped with two 2° wedged 
polyethylene (PE) windows to minimize channeling fringes. The 
results from these samples were used as a supplemental that, in 
combination with the waveguide data, allows determination of 
the clump size of our particles. 

Dust diluted in IR-transparent matrix materials, such as 
PE and KBr, are used for shorter (mid-IR) wavelengths. For 
these samples, the dust is mixed into the substrate powders 
and then pressed (KBr samples) or melted (PE samples) into 
a disk. The amount of dust chosen to be mixed into the 
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substrate is small enough to allow it to be optically thin for 
transmission measurements. Additionally, the preference is for 
the spectra from various wavelength regimes to overlap when 
comparing the optical properties of the dust measurements. 
The dust concentrations are varied for each disk to provide 
moderate (20%-80%) transmission across the wavelength range 
of interest. Reference disks with no dust sample material were 
used to calibrate the measured dust spectra. 

Moderate resolution (4 cm -1 ) transmission spectra for the 
samples are measured in the mid- and far-IR (10-4000 cm -1 or 
1000-2.5 /rm) using Bruker IFS 125HR and 113v spectrome- 
ters. Each spectrum was produced from 64 scans with the FTS 
to allow for a statistical spectrometer error of less than 1%. 
For temperature-dependent studies (i.e., the aluminum holders 
packed with dust and the dust diluted in PE), the sample is 
mounted in an Oxford Optistat continuous-flow liquid helium- 
cooled cryostat equipped with an Oxford ITC503 temperature 
controller in order to measure the sample transmission at de- 
sired temperatures from 5 to 300 K. Alignment repeatability 
problems within the cryostat for smaller-diameter KBr disks, as 
well as the hydroscopic qualities of KBr, limited our capabilities 
to perform cryogenic measurements in the mid-IR. 

2.4. FTS Reflection Measurements 

To supplement the transmission measurements, diffuse reflec- 
tivity in the mid- and far-IR is measured using a novel reflec- 
tometer apparatus, operable at temperatures below 100 K. The 
system is equipped with an integrating sphere, a Si bolometer 
detector element, and a three-position sample wheel thermally 
isolated from the helium bath. The sample wheel includes po- 
sitions for the sample, a polished aluminum disk for reference, 
and a black sample calibrator. The dust sample is deposited on 
an aluminum disk by first diluting the dust in acetone and then 
depositing the solution onto the A1 sample disk. After the ace- 
tone evaporates, a relatively homogenous dust coating remains 
on the A1 disk. This homogeneity over the size of the beam is 
such that the spectra are reproducible for several samples. The 
operation of the reflectometer is detailed in Rinehart et al. (20 1 1 ) 
and the black sample calibrator in Quijada et al. (201 1 ). The re- 
flectometer can be used to measure reflectivity at 100-650 cm -1 
(100-15.4 fi m) using the Bruker IFS 1 13v. 

Mid-IR (400-1200 cm -1 , 25-2.5 fim) reflection data is taken 
using a gold-plated integrating sphere accessory available for 
the Bruker IFS 125HR. A focused beam is directed toward the 
sample with a 13° angle of incidence and is read by a DLaTGS 
detector equipped with a KBr window. As this accessory is not 
capable of being used cryogenically, only room temperature 
reflection data were measured. 

3. MODELING AND ANALYSIS 

The real and imaginary components of the optical function 
(the refractive index, n, and absorption coefficient, k respec- 
tively) are frequency-dependent and determine how the radi- 
ation is transmitted and reflected as it propagates through a 
sample. The analysis of the laboratory spectra employs a set 
of mathematical models for the dielectric slab samples that are 
applied according to the degree of optical coherence within the 
sample under the assumption that the mixture is homogenous 
and the sample geometry has planar faces which are parallel. A 
general one-layer slab model (Bohren & Huffman 1998) is ap- 
plied to compute the resulting optical response for plane-parallel 
interfaces in the limit the particle size small compared to a 


wavelength and the coherent reflection and transmission oc- 
cur. In the limit the wavelength is large compared to the grain 
(particle) size, the channel spectra fringes fade due to increas- 
ing decoherence and subsequent scattering of the light out of 
the FTS beam. In this limit, the transmission and reflection re- 
sponses can be characterized by treating the optical parameters 
as a power law (Halpem et al. 1986). More generally, the scat- 
tered light in a plane-parallel sample of homogenous material 
are partially coherent and can be treated by the approach de- 
scribed in Grossman & McDonald ( 1 995). 

3.1. Waveguide Measurements 

At microwave frequencies (~ 1 cm -1 ), we model the response 
of the waveguide sample holder as a multi-section transmission 
line via the “ABCD”-matrix approach (Pozar 2004); this tech- 
nique is also known as the chain, transmission line, or char- 
acteristic matrix in the literature (Goldsmith et al. 1998; Yeh 
1988). Physically, this matrix formulation links the propagation 
of the input and output electric and magnetic fields in a layered 
medium with plane parallel boundaries. This computational ap- 
proach enables the perturbative influence of the windows to be 
incorporated into the modeled Fabry-Perot resonator response 
by forming the product of the “ABCD”-matrix for each element 
of the sample holder. It is useful to recall that the propagation 
constant, y, for a waveguide section is given by, 



where a> is the angular frequency of the incident radiation, c is 
the speed of light in the vacuum, and the cutoff wavelength X c is 
two times the waveguide broadwall length for the dominant 
mode in a rectangular waveguide. The complex dielectric 
permittivity and magnetic permeability of the material filling 
each section of the guide are parameterized as 


e = e*e 0 - (e' r + ie")e 0 . 

(2) 

v- = /^>o = 04 + //4Vo> 

(3) 


where to 8.854 x KT 12 (F m _1 ) and /z 0 = 4n x 1CT 7 
(H m '). The analytical expressions in Wollack et al. (2008) 
provide a useful limiting form where windows are not present 
or their thickness relative to the guide wavelength tend toward 
zero and can be neglected. Using Equations (2) and (3), we can 
determine the refractive index as 

n* = n + ik = y/efji*, (4) 

where e* is the complex relative permittivity and ti* the complex 
relative permeability. From this we find 


_ I J(K) 2 + (£'r) 2 + K 

V 2 


and 


k = 


_ I y/(g;) 2 +«) 2 -e; 
2 


(5) 


( 6 ) 
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3.2. FTS Transmission Measurements 

For the matrix sample measurements (KBr or PE), sophis- 
ticated models are required (Sihvola 1999). The dust particles 
are embedded in a background material with a small volume 
filling fraction, allowing the samples to remain optically thin at 
shorter wavelengths. The mixture is considered dilute, and the 
Maxwell-Gamett formula can be used to describe the dielectric 
properties (Garnett 1904). The effective dielectric function, % f , 
is expressed as a function of the volume filling fraction, /(typi- 
cally 10 3 ) and the dielectric functions for the background and 
inclusion, fiback and £i nc . respectively: 


Sets — ®back + 3 f£ back 


£inc ^back 

£inc + 2£“back f{£inc “ ^back) 


(7) 


For KBr and PE, fback can be approximated; we used the 
value found in Palik (1985). £i nc can be expressed as a set of 
Lorentzian oscillators: 


M 


Sine = (n + ik) 2 = 


inc.oo 


■5> 

j = i 


to 


p-j 


a>lj — to 2 — icoSj 


( 8 ) 


Each oscillator is defined by the variables o> 0 (the resonant fre- 
quency), t o p (the plasma frequency), bj (the oscillator strength), 
and S (the full-width at half-maximum). This enables the diluted 
materials to be parameterized and solved from the data obtained. 
%r can be inserted into the equation for transmission (Bohren & 
Huffman 1998), utilizing the relationship between the complex 
effective dielectric function and the complex effective index of 
refraction, n c( t + for the mixture: 


where 


and 


(1 - ftefr) 2 exp(-g eff h) 

1 - fl| ff exp(— 2aefrfi) ’ 


*cff = 


(rteff 1 ) 4" A'^|y 
( r hl! + l) 2 + k\ if ' 


= 47T vA'eff, 


(9) 

GO) 

(ID 


where a e n is the effective absorption coefficient. There are 
models that may improve the quality of the fit of broader features 
(e.g., the feature caused by water in the range 2600-4000 cm - 1 ), 
by introducing an adjustable Gaussian (Brenden & Bormann 
1992; Cataldo et al. 2012). The current models are sufficient for 
use with these data, as the remaining residuals after accounting 
for water in the sample are less than 5% (compared to 10%-15% 
residuals in fitting the dust parameters). However, the improved 
models are in development for use in interpreting higher signal- 
to-noise data obtained in the future. 


3.3. FTS Reflection Measurements 

The analysis of the reflection data uses a modified version of 
Beer’s law, which accounts for the scattering component of the 
absorption coefficient (Egan et al. 1973): 


and 


(1 - /? 2 )e -r 
1 - R 2 e ~ 21 ’ 


( 12 ) 

(13) 
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where 

0 = J~’ (14) 

V a + s ' 

and 

r = 2y/a{a + s)h, ( 15 ) 

where s represents the total scattering coefficient, h the thickness 
of the sample, a is the absorption coefficient, and t is equal 
to oferrh. These equations are valid for the case of the diffuse 
transmitted component of the radiation. 

The mathematical models for the matrix-transmission mea- 
surements were implemented in MATLAB through a least- 
squares nonlinear fit of the transmission equations to the lab- 
oratory data. The solver is a Levenberg-Marquardt algorithm 
with finite-difference computation of the Jacobian matrix. Sim- 
ulations were run forcing the reduced x 2 to be close to unity, 
in order to estimate the systematic error affecting the measure- 
ments. This was found to be ~2%, which is within our expected 
accuracy range. The parameters are recovered to an accuracy of 
10 -6 from simulated transmittance spectra (Cataldo 2010). 

Due to the fluffy, open aggregate nature of these grains, the 
density is not definitively known, however, our models require 
a value for parameterization of the optical functions. In order 
to test how the density will affect the outputted n and k values, 
we ran our model for the 52 mg iron silicate in a 1000 mg PE 
matrix for 4 different densities: l,2,3,and4g cm -3 . Then value 
was the same for all four densities at 435 cm -1 , and the most 
extreme case of differing values occurred at 40 cm -1 , where 
the p = 4 g cm -3 n value was 1.24 times stronger than the 
P = 1 g cm -3 n value (with p = 3 g cm -3 and p = 2 g cm -3 
in between). The k value showed a greater discrepancy in the 
outputted values between different densities. Again, the extreme 
value occurred at 40 cm' 1 , with the k value for the p = 4 g cm -3 
being 4.4 ( k = 0.44) times stronger than the p = 1 g cm -3 
value ( k = 0.10). Studies of interplanetary dust particles have 
indicated densities ranging between 0.3 and 6.2 g cm -3 , with 
iron-containing particles near 3.5 gem -3 (Love et al. 1994). 
Due to the iron-rich nature of our samples, we have estimated a 
density of 3.8 g cm -3 , which is slightly lower than the average 
density of ferrosilite (3.95) or fayalite, and iron-rich olivine, 
with a density of ~4.4 g cm -3 . 

4. RESULTS 

4. 1. Waveguide Measurements 

The spectra and modeled fit for the waveguide data are 
shown in Figure 2. For the purposes of fitting the observed 
data, the measured waveguide length, broadwall, and window 
parameters were fixed and the model’s complex permittivity and 
permeability were varied. The deviation between the modeled 
and observed reflection and transmission data sets are jointly 
minimized during the fit. The details of the window material 
between the waveguide flanges are presently untreated in the 
model. This influence can be seen as a resonance dip in the 
average of SI 1 and S22 near 28 GHz, however, given the joint 
treatment of the scattering parameters, the extracted parameters 
are verified not to be sensitive to this portion of the data at 
the reported precision. For the iron silicate sample investigated 
here this is not a limiting concern. Should this systematic effect 
become an issue in the future with different types of samples, 
it can be moved out of the sample holder’s effective pass-band 
either by reducing the window’s electrical length by a factor 
of approximately two or by eliminating it using a photonic 
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WR28.0 Sampleholderr Iron Silicate 



Figure 2* Using a waveguide resonator packed with an iron silicate sample, we derive the complex dielectric function at millimeter wavelengths. Shown are the 
measurements for loosely packed material as well as our model lit, indicating a clean lit with the data (error within a few percent). The model, combined with 
knowledge of the lilting fraction of the material allow us to calculate n - 2.7 and k = 0,20 for the bulk iron silicate material from 22 to 40 GHy. (averaged at M cm -1 ). 
(A color version of this figure is available in the online journal) 


waveguide choke flange on the VNA reference port interface 
(Wollack el ah 2010). 

From the observed spectral response of the iron silicate dust 
in the waveguide sample holder, £*(dust) = 1.067 + / 0.01 1 
(/i e ft — 1.034 and = 0,053) and /f*(dust) = 1.0 from 
22 to 40 GHz (—1 cm" 1 ). Taking p = 3,8 g cm -3 for the 
bulk density and the sample’s estimated 1 .8% volume filling 
fraction, £*(bulk) = 7.2 + / 1.1 (nj [lc ~ 2,7 and k inc — 0.20) for 
a hypothetical bulk dielectric sample of this material. 

Utilizing the f*(bulk) calculated from the waveguide, we 
utilized iVlie theory with transmission spectrum from the 2 mm 
and 4 mm aluminum holders in the far-lR to estimate the dump 
size of the particles, shown in Figure 3, A Mie theory estimate 
of the seatting size from the FTS in section gives a calculated 
clump radius of 19,4 /mi at room temperature. This estimate is 
within the assumptions used in fitting waveguide data. Note that 
as the temperature decreases, the clump size increases from a 
clump size of 19,4 p m at room temperature to a clump size of 
24,0 /nn below 1 00 K. The maximum calculated error for the fit 
of the room temperature data was 5%, and tins value increased 
as the temperature decreased (to a value of 1 2% for 50 K), The 
average calculated error for each fit was -~2%. 

4.2. FTS Transmission Measurements 

Figure 4 displays the transmission spectrum for FeSiO 
condensate diluted with a substrate of melted PE for the 
40-400 cm -1 range (250-25 pm) at room temperature. For this 
sample, 52 mg of condensate were mixed into 1000 mg of PE 
pellets and then melted into a solid disk; the beam measuring the 
sample passes through the center of the disk where the thickness 
is ^1.25 mm and the computed filling fraction is 1.23 x 10“\ 
Figure 5 displays the calculated optical parameters, n and A r , 
for the dust diluted in PE, The n across the spectrum varies 


significantly with frequency (between ~2,6 at lower frequencies 
and M.7 at 400 cm" ). 

Data was also collected tor the FeSiO condensate diluted 
with a substrate of melted PE at 300, 100, 50, and 5 K 
for 100-650 cm" 1 (100-15,4 /Am). For this sample, 15 mg 
of condensate were mixed into 1000 mg of PE pellets and 
then melted into a solid disk; the beam measuring the sample 
passes through the center of the disk where the thickness is 
M .78 mm and the computed filling fraction is 3.58 x J0“ 3 and 
the spectrum is shown in Figure 6. The spec mini is dominated by 
a strong absorption feature, which has a minimum at 464 cm -1 
(21.5 /nn), likely resulting from a Si-O rocking mode (Boudei 
et al. 2005). The spectral region dominated by the absorption 
band shows no change in transmission with temperature and 
only a small amount of deviation at lower frequencies (e.g.* 
— 5 % difference at 100 cm -1 ); the majority of this deviation 
occurs between 300 K and 100 K, with little additional change 
below 1 00 K, Variation of n and k with temperature are similarly 
small, as shown in Figure 7, Differences in the value of n are only 
noted at the low frequency end of the range, while the values 
of A show somewhat larger variation between room temperature 
and cryogenic temperatures. 

Mid-IR data for the dust in a KBr mixture was obtained 
using a dust- KBr pellet. Figure 8 displays the transmission 
spectrum for 0.59 mg of FeSiO condensate pressed into a 
pellet of 500 mg KBr from 400 to 4000 cm" 1 (25-2.5 pm). 
The spectrometer beam passes through the center of the disk 
where the thickness is 1 ,55 mm and the computed volume filling 
fraction is 4.05 x 1 0 -4 . Data for this sample were collected only 
at room temperature as described in the Experimental Methods 
section. The spectrum includes the 464 cm -1 (20.6 pm) feature 
from a Si-O rocking mode thai was also observed in the dust-PF 
sample. Additionally, the 800 cm -1 (12.5 gm) feature from the 
Si-0 bending mode and the 1030 cm" 1 (9,7 pm) feature from 
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Wavelength [pm] 



Figure 3. The inmsirasskn spectrum for the iron silicate Uust in 2 mm and 4 mm aluminum holders at room temperature (black), as well as the fit from our model 
(red). Using Mie Theory to fit these data allows us to determine a room-temperature particle clump %\7£ of 19.4 /ini (with an average error of 1%), 

(A color version of this figure is available in the online journal.} 


Wavelength [pm] 



Wavenumbers (cm' 1 ] 

Figure 4, The transmission spectrum for 52 mg of iron silicate condensate diluted with a HHM) mg substrate of melted polyethylene for the 40-400 cirT range at 
room temperature (black), as well as the lit from our model (blue) and the calculated residual of the fit (red)* 

(A color version of this Jigurc is available in the online journal .) 


the Si-O stretch band are shown. These feature positions arc 
consistent with those in other studies ( Bonder et al. 2005). 

Both the iron silicate sample and the KBr matrix material 
are hygroscopic; therefore, the models used for analysis were 
adapted to account for the presence of water in the spectrum. 
Water trapped in the matrix is identified by the O-H stretch- 
ing mode at 3280 cm -1 (3,05 /an), the O-H bend mode at 
1660 cirT 1 (6.02 /im), and the libration mode at 760 cm -1 
f 13.2 /tm). This complicates the Lorentzian model and the fit- 
ting process. In order to ensure we utilized the proper filling 
fraction, the optical functions for water were compared to pi e- 
vious studies (Hale & Querry 1973) and were found to be in 
excellent agreement. Taking the effects of water into account 
in the fitting process provides an improvement in the n and k 
analysis compared to previous attempts to model this type of 


data (K ravels et ah 2005). As shown in Figure 8. the normalized 
residual is well below 10% for the majority of the spectra. How- 
ever, the normalized residual for the fit (defined as the absolute 
difference between the experimental and fit spectra) is higher on 
average than other data analyzed within this paper. One spike 
in the residual occurs between J238 and 1290 cm' 1 , reaching 
a peak normalized residual at 15%, This spike in the residual 
value results from the relatively low amount of transmission, 
followed by a drastic increase in the slope of the curve. When 
the slope is large, the variations between the fit and the exper- 
imental data is exasperated. This, combined with the difficulty 
in fitting transmission data below 0,2. results in the higher error. 

Figure 9 displays the calculated optical parameters, n and 
k, for the iron silicate dust diluted in KBr. The n across the 
spectrum varies between 2.6 and 0.60, trending to l .49 at higher 
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Figure 5, From ihe transmission spectrum shown in Figure 4, the optical parameters, n (blue) and it (red), can be calculated for the iron silicate inclusions, as outlined 
in the modeling and analysis section from 40 to 400 em~ ] for iron silicate dust diluted in polyethylene when modeled as a multi -constrained function. 

(A color version of this figure is available in the online journal) 


Wavelength [pm] 

16 20 25 30 40 50 65 100 



Figure 6. The transmission spectra for 15 mg of iron silicate condensate diluted with a 1000 mg substrate of melted polyethylene from ihe 100-650 cm 1 range at 
room temperature (red). 100 K (green), 50 K (yellow), and 5 K (blue). 

(A color version of this figure is available in the online journal.) 


frequencies. The k values exhibit two distinct peaks (1.2 at 
448 cm -1 and 1.7 at 1062 cm” 1 ); at higher frequencies, k lends 
toward zero. 

4.3. FTS Reflection Measurements 

For ihe samples described in Section 2,4, supplemental 
diffuse reflection data was collected at 10, 50, and 100 K for 
the wavelength range of 400-700 cm" 1 (25-143 /mi). As ihe 
MX4S5 mm sample was optically thick, there was no significant 
diffuse reflection at any of those temperatures. For the optical 
functions, the refractive index was calculated to be n — 2.2 and 
the absorption coefficient was calculated to be k — 0.5 near 
200 cm -1 and 0.9 near 400 cm -1 . 

Additionally, the refection data obtained in the mid-IR, in the 
range of 400-1200 cm” 1 and at room temperature, did not show 


significant diffuse reflection. For the optical parameters, the 
refractive index was calculated to be n = 1,4 and the absorption 
coefficient was calculated to be k ~ 0. 1 1-0,15. Because the 
measured reflection from the samples was very low, model fits 
necessarily produced large residuals and had large associated 
uncertainties. Therefore, both ihe far-IR and mid-lR reflection 
results were only used in comparison with iransmission data and 
for verification purposes. 

4.4. Comparison between Sample Sets 

In order to understand the data over the entire wave- 
length range, comparisons of ihe optical properties between 
different sample sets arc made and shown in Table 2 and 
Figure 10. Regions of overlap between samples typically 
showed more variation, due to limiting factors within our 
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Wavenumbers [cm -1 ] 

Figure 7, From the transmission spectra shown in Figure 6, the optical parameters, n and k, can he calculated for the iron silicate inclusions, as outlined in the modeling 
and analysis section from KX) to 650 cm” 1 for iron silicate dust, diluted in polyethylene. The refractive index (n) for the iron silicate inclusions in the samples over 
various temperatures overlaps nicely* while the absorption coefficient (it) shows only a small v ariation with temperature. Note, due to a lack of noticeable change with 
temperature, the 100 K and 50 K spectra are difficult to sec and lie under the 5 K spectra, 

( A color version of this figure i s available in Lhe online journal*) 


Wavelength [him] 


25 3 4 5 10 20 



Figure 8* The transmission spectrum for 0.59 me of iron silicate condensate presM;d with a 5(H) mg pellet of KHr over the 400-4000 cm” 1 range at room temperature 
(black). Water trapped in the matrix complicated the Loren Loan model and the lit; taking this into account, the model lit shown (blue), as well as the residual in the fit 
(red). 

{A color version of this figure is available in the online journal) 


Table 2 

n and k Values for liach Sample Measured 


Sample Type 

Wavenumber Range 

Range of n Values 

Range of A Values 

FeSiO in waveguide 

M cm -1 

2*7 

0.20 

52 rng FeSiO 1000 mg PF 

40-400 cm” 1 

1. 7-2.6 

0,494)93 

15 mg FeSiCN-1000 mg PF 

100-650 cm" 1 

1 .3-2.3 

0*25-0*87 

0.59 mg leSiO+50G mg Klir 

400-4000 cm” 1 

0*56-2.5 

<0*01-1.6 

Reflection: FcSiO+aeetone 3 

100-650 cm 

2.2 

0*02-1*5 

R e flee tion : I ; eSiO+ace tonc JI 

400-1 200 cm" 1 

1,4 

0.11-0*15 


Notes* 1 Used for verification purposes only. See Section 4,3. 
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Figure 9. From the transmission spectrum shown in Figure X. the optical parameters, n (blue) and k (red), c;m bo calculated for the iron silicate inclusions, as outlined 
in the modeling and analysis section, for 400^4000 cm - 1 for iron silicate dust, pressed within a Kl Jr pellet. when modeled as a mu In -con strained function. 

(A color version of this figure is available in the online journal,) 



Figure 10. The calculated optical parameters, n (blue) and k (red j. for the iron silicate inclusions from 1 to 4000 cm -1 . Multiple sample types and mixtures were used 
in order to cover the wavelength region accurately. Regions of overlap between samples typically showed more variation, due to limiting factors within our fitting 
procedure for the extreme upper and lower wavelength regions of each sample. 

(A color version of this figure is available in the online journal.) 

lining procedure for the extreme upper and lower wavelength 
regions of each sample (typically also shown as an increase in 
error for the fit). However, as shown in Figure K), the general 
trend of the optical parameters for ihe iron silicate inclusions is 
visible. 

For the millimeter wavelengths (—1 cm" 1 ), ihe refractive 
index, n, was calculated to be 2.6, Although the two samples do 
not have an overlap in wavelength regions, the dust- PE sample 
shows an increase in n approaching the millimeter and could 
be used as an indicator for our sample converging toward that 
value. The 52 mg FeSiO + 1000 mg PE sample was measured 
lo values further into the far-IR (to 40cm -1 ) for verification 
and continued to show an increase in the refractive index 
values that could further lead to justification of these values. 

The 52 mg FeSiO + 1(300 mg PE sample overlaps well with 


the 15 mg FeSiO-r 1000 mg PE sample in the 150-400 cm -1 
range, as shown in Figures 5 and 7. The n values for the 15 mg 
FeSiO + 1 000 mg PE and 0.59 mg FeSiO + 500 mg KB r samples 
overlap very well in refractive index, as shown in Figures 5 and 9. 
In the mid-IR* the refractive index is calculated to be 1 .4, which 
aligns well with the dust-PE measurement trend toward 1 .4 in 
the tnid-IR. In relation to the mid-IR data, the refractive index 
again aligns well with the 0.59 mg FeSiO + 500 mg KBr sample, 
with the n value trending toward 1.5 (only a —5% difference 
from the mid-IR reflection data). The n from the reflection data 
in the far-IR was determined to he 2.2; this is consistent with the 
range of values found from the dust-PE samples* which range 
between 1.7 and 2.6 in the far- JR. 

The absorption coefficient* k, varies more between samples in 
comparison to the values for the refractive index. In the far-IR 
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(to 10 cm” 1 ), the 52 mg FeSiO+ 1000 mg PE sample is slightly 
higher than those for the waveguide measurements, however, 
we would expect a downward trend to occur between 40 and 
1 cm -1 . The alignment of the 52 mg FeSiO +1000 mg PE 
sample and the 15 mg FeSiO+ 100 mg PE sample is good in 
the lower wavenumber region, however, the k values increase 
slightly more in the 15 mg sample above 300 cm -1 , This is due 
to the transmission of the 52 mg FeSiO sample being low enough 
in this portion of the spectrum. For the 15 mg FeSiO + 1000 mg 
PE sample and the 0.59 mg FeSiO + 500 mg KBr sample, the 
alignment is good, as shown in Figures 7 and 9. Additionally, 
for the far-IR measurements, the reflection measurement shows 
a larger distribution of k values; between 200-450 cm” 1 , the 
data aligns within a factor of 15% (with k values around 0.5 
near 200 cm -1 and 0.9 near 450 cm” 1 ). The results from all of 
these measurements are consistent, taking into account the large 
uncertainties associated with the reflection data. 

5. DISCUSSION AND IMPLICATIONS 

By using the calculated optical functions shown here, model- 
ing of thermal dust emission spectra can be accurately performed 
within the wavelength range covered by our results. Due to the 
unique nature of our iron silicate samples, our experimental 
data can be used for comparison with observations from current 
and future missions, including Herschel , Spitzer , SOFIA, JWST, 
and ALMA. However, in comparing our results with astrophys- 
ical environments, we need to fully understand how silicates 
behave within those environments. In particular, while results 
from Nuth et al. (2002) showed that pure, amorphous iron sil- 
icate grains condense from a mixed Fe-Mg-silicate vapor, the 
magnesium silicate grains anneal at much lower temperatures 
than the iron silicates. Therefore, iron silicate grains will remain 
amorphous for longer periods of time than the magnesium sili- 
cate grains. While the magnesium silicate grains can be utilized 
for calculation of the mass-loss rate, the slower processing of 
iron silicate grains makes them unsuitable for tracing the history 
of dust outflows. However, iron silicates are useful in determin- 
ing the abundance and enrichment of dust grains, and is clearly 
needed for mixtures of Fe-Mg grains and for understanding the 
mass-loss history. 

6. SUMMARY AND FUTURE WORK 

The OPASI-T program utilizes multiple instruments to pro- 
vide spectral data over a wide range of temperatures and wave- 
lengths of dust grain analogs. Experimental methods include 
waveguide, transmission, and reflection/scattering measure- 
ments. From this data we can determine the optical functions, n 
and k. The analysis of the laboratory data for each sample type 
is based upon different mathematical models, which are applied 
to each data set according to their degree of coherence. In the 
work presented here, iron silicate dust grain analogs (in several 
sample types and at temperatures ranging from 5 to 300 K), 
were studied across the IR and millimeter portion of the spec- 
trum (from 2.5 to 1 0,000 /zm or 4000 to 1 cm" 1 ), and the optical 
properties, n and k , were calculated. Fits of the models to the 
iron silicate data were very good, particularly for the waveg- 
uide data in the millimeter and the melted dust-PE disks in the 
far- and mid-IR. Water complicated the fit for the dust-KBr 
disk in the mid-IR, which lead to a higher error, but still us- 
able optical property calculations. The far-IR transmission data 
using the aluminum holders allows determination of the clump 
size, and the reflection data in the mid- and far-IR was used for 


verification purposes with the data. Experimental measurements 
have been taken for magnesium silicate samples, data analysis 
to obtain optical functions has been completed, and the publi- 
cation is currently in preparation. Additionally, annealed mag- 
nesium silicate samples have been produced and have almost 
completed the process of spectral collection and data analysis. 
Comparisons will be made to the iron silicate data, as well as to 
previous works in literature (e.g.. Speck et al. 201 1). 
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